Unified description of Fermi and non— Fermi liquid behavior in a conserving slave 
boson approximation for strongly correlated impurity models 
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We show that the presence of Fermi or non-Fermi liquid behavior in the SU(N)xSU(M) An- 
derson impurity models may be read off the infrared threshold exponents governing the spinon 
and holon dynamics in a slave boson representation of these models. We construct a conserving 
T-matrix approximation which recovers the exact exponents with good numerical accuracy. Our 
approximation includes both coherent spin flip scattering and charge fluctuation processes. For the 
single-channel case the tendency to form bound states drastically modifies the low energy behavior. 
For the multi-channel case in the Kondo limit the bound state contributions are unimportant. 
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Impurity models with internal degrees of freedom and 
strong local correlations coupled to a fermionic bath have 
been of considerable interest recently The proto- 

type is the Anderson impurity model, involving a local- 
ized electron level (called <i-level in the following) hy- 
bridizing with one or several conduction bands ||. The 
strong Coulomb repulsion U (U — > oo) between electrons 
in the localized state effectively restricts the d-level occu- 
pancy to rid < 1- The ensuing projection of Hilbert space 
onto the physical subspace without multiple occupancy 
is a problem of fundamental importance in the theory 
of strongly correlated Fermi systems in general. As a 
consequence, the Anderson model displays many of the 
salient features of strongly correlated systems, including 
the formation of local magnetic moments and a compe- 
tition between non-FL behavior caused by an incipient 
orthogonality catastrophy, to which the system scales ini- 
tially, and a FL fixed point, which is realized at energies 
below a characteristic scale, the Kondo temperature Tk, 
if the local moment can be completely screened by the 
conduction electron spin system. This model can, there- 
fore, serve as a test case for the regime of strong correla- 
tions and at the same time for developing new methods 
which may later be applied to lattice problems as well. 

In terms of pseudofermion and slave boson operators 
U, b m (a = 1,...,N, m = 1,...,M) (|f] the M- 
channel Anderson model is defined by the Hamiltonian 

H = H + E d flU + V ]T (4 blU + h.c), (1) 



where H a = V\- (er 
duction band energy. 
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creates a conduction elec- 
tron in band m with spin projection a and momentum 
k, Ed denotes the energy of the AT-fold degenerate lo- 
cal d-level at position R = 0. V is the hybridization 
matrix element and fi the chemical potential. A phys- 
ical electron in the local level is created by the elec- 
tron operator dy = ^2 m f„b m , where the condition of 
no double occupancy is effected by the local operator 



constraint Q = J2 a faf<y + J2 m b ln b m = 1- The effec- 
tive coupling constant in the constrained Hilbert space is 
given by r = irV 2 N(0), with N(0) the conduction elec- 
tron density of states at the Fermi level. For M > 1 the 
Hamiltonian Eq. ([[]) does not correspond directly to a 
physical system, since it involves the existence of several 
empty orbital states b^vac), while there is usually only 
one in a physical system. Introducing these states is a 
mathematical convenience which allows one to derive the 
SU(N)xSU(M) Coqblin-Schrieffer model from Eq. © in 
the Kondo limit. The latter has been studied extensively 
by the Bethe ansatz method M], conformal field theory 
(CFT) pi and self-consistent slave boson theory jq] . 

In this Letter we focus on the auxiliary particle Green's 
functions G /<r (n -r 2 ) = -(T{U(n)fifo)}), G b m (n - 
t 2 ) = — (T{6 m (ri)6j 7i (T 2 )}). The angular brackets denote 
the statistical average in the grand canonical ensemble, 
((. . .)) = tr[(. . .)exp(-/3(#-AQ))]/tr[exp(-/?(Zf-AQ))]. 
The exact projection onto the subspace Q = 1 is achieved 
by differentiating with respect to the fugacity exp(— (3X) 
and taking the limit A — ► oo 11 1. This procedure 
ensures that the projected propagators obey Wick's the- 
orem, and self-energies 'Sf t b,c(i^>n) m ay be defined by 
Gf t b,c(iu n ) = ([G/ ibiC (w n )] _1 - S/ )6)C (iw n )) , where 
iw n ) = (icj n -E d - A)~\ G° bm (iuj n ) = (ioj n - A)~\ 

1 The d electron 



and G° cma (iu)n) = £s0 
Green's function may be expressed in terms of E c as 
discussed in Ref. Q. The projected spectral functions 
A x {u>) = lmG x (oj — iO) exhibit divergent threshold be- 
havior at ijj = with a proper choice of the zero of the 



auxiliary particle energy 111 



A x (lo) (xuj Qa! , x = f,b. 



For the single-channel model with spin degeneracy N, 
which is known to have a FL ground state, the exact ex- 
ponents a.f t b have been determined ]12 j by Wilson's nu- 
merical renormalization group (NRGj approach for the 
cases N = 1,2, M = 1. They may also be deduced for 
arbitrary N by the following argument (i) In the 

spin screened FL state (u>, T < Tk) the impurity is seen 
by the conduction electrons as a pure potential scatter- 
ing center, (ii) The infrared (IR) threshold behavior of 
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Gf.b is then entirely due to the orthogonality catastro- 
phe of the overlap of the Fermi sea without the impurity 
level and the fully interacting conduction electron sea af- 
fected by the potential scattering phase shifts 8 a . The 
corresponding exponent is given by a = 1 — ^2 (T (S (T /tt) 2 . 
(iii) The phase shifts follow from the Friedel sum rule 
An a = 5c /n, where An a is the change in the number of 
conduction electrons at the impurity caused by the in- 
teraction with the impurity, (iv) For the boson spectral 
function the initial state is the empty impurity, which 
for each spin species fills up with An a = n^/N conduc- 
tion electrons in the final state, until the correct impurity 
level occupation rid is reached. It follows that 

a b = l- n 2 d /N. (2) 

For the spectral function of fermions with spin a the ini- 
tial state is defined by a full impurity level with spin a 
with the remaining N — 1 impurity levels empty. The 
corresponding change of conduction electron number in 
the final state with occupation nd is An a — nd/N — 1 
and An a i — n d /N , a 1 ^ rr, and hence 

at = (2n d - n 2 )/N. (3) 

We emphasize that the expressions Eqs. (||),(||) for the 
exponents have been confirmed using the Bethe ansatz 
solution and boundary CFT |i~4]]. In the Kondo limit 
(n d = 1), cif = 1/N, in disagreement with a result de- 
rived from a self-consistent parquet analysis |L5|. Note 
that complete spin screening is crucial for this derivation 
of the exponents in terms of scattering phase shifts to be 
applicable: For the multi-channel model (M > 1), which 
exhibits a non-FL ground state, the exponents in the 
Kondo limit are known from CFT [§ to be a f = M/(N+ 
M), ab = N I (N + M), while the above argument would 
yield the wrong result a/ = 1 — M + (2nd — n d /M)/N, 
at, = 1 — n 2 d /(NM). Therefore, the IR threshold expo- 
nents of the auxiliary particles are indicators, i. e. a nec- 
essary and sufficient condition, for FL or non-FL behav- 
ior, respectively. 

There is evidence that a/ has also observable rel- 
evance in that it governs the physical electron spectral 

function at intermediate frequencies ui ~ Tk- 

We now turn to an approximation scheme |l7j which 
is capable of recovering the above (exact) IR dynamics. 
As a minimal requirement, the constraint Q = 1 has to 
be fulfilled in any approximate theory. The constraint 
is closely related to the invariance of the system un- 
der a simultaneous local (in time) gauge transformation 
/ ct (t) ->■ exp(i9(T))/ (T (r), 6 m (r) ->■ exp(i6(r))6 m (r). 
The Lagrange multiplier A assumes the role of a lo- 
cal gauge field and transforms as A — > A — id® /dr. 
Any approximation scheme respecting the gauge sym- 
metry will preserve the charge Q in time. We shall 
call approximations of this type conserving. Symme- 
try conserving approximations of the self energies E;>./, c 
and the irreducible vertices T xy may be generated by 
functional differentiation from a functional $ of closed 



skeleton diagrams as ^^(ri — t^) = 5$>/8G x (ri — T2), 
r xy = 5 2 $/(6GJGy), x,y = f,b,c. 

We will be interested in the limit of weak hybridiza- 
tion, such that the dimensionless parameter VN(0) <C 1. 
Thus, let us first discuss the lowest order approxima- 
tion, which is of second order in V. The conserving 
approximation scheme requires the self-energies to be 
determined self-consistently, which amounts to an in- 
finite resummation of perturbation theory even if only 
the second order skeleton diagram for <E> is kept. The 
resulting scheme is known as the "Non-crossing approx- 
imation" (NCA) p)-|0|. It should be a qualitatively 
correct approximation, provided the perturbation series 
for <f> converges, i.e. if there are no additional collective 
effects causing singularities. The NCA leads to very 
good results in the absence (i.e. in the multi-channel 
case) or sufficiently far away from a FL fixed point: 
A comparison of NCA results for the auxiliary particle 
and ci-electron spectral functions Af, Ab, Ad and ex- 
act results obtained for the single-channel case using the 
NRG method shows that (i) the NCA auxiliary par- 
ticle spectral functions are even quantitatively correct 
at energies lo around and above the Kondo scale Tk = 
[MT /ity M l N '>eyLv[-it\E d \/(NT)], but (ii) their low en- 
ergy behavior (lo Tk) is incorrect. The latter appears 
to be due to a lack of vertex corrections. Within NCA the 
exponents of the above mentioned threshold power laws 
may be determined analytically as aJ CA = Mj (N + M), 
a NCA = N /( N + M ) For the case M = 1, 

these values disagree strongly with the exact results dis- 
cussed above. In the multi-channel case (M > 1), on 
the other hand, the NCA exponents agree with the ex- 
ponents found for the fundamental fields and their corre- 
lation functions in CFT in the Kondo limit. This sug- 
gests that the NCA describes the low energy properties 
correctly in the non-FL regime of the SU(N) x SU(M) 
Anderson model for nd = 1 and that the generic behavior 
of the model is that of a non-FL. 

It may be shown by power counting arguments that 
there are no corrections to the NCA exponents in any 
finite order of perturbation theory [||. However, addi- 
tional collective effects, e.g. the formation of the Kondo 
singlet state, lead to FL behavior. Thus, it is natural to 
search for singularities in the pseudofcrmion-conduction 
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FIG. 1. a) Diagrammatic representation of the conduction 
electron-pseudofermion T-matrix T^h b) Pseudofermion 
and c) slave boson self-energies E/ CT , £;,. The terms corre- 
sponding to T^ 6 ' are obtained by interchanging / <-» b. 
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electron scattering channel. In particular, we consider 
the class of diagrams which, at any given order of V 2 , 
represents processes with the maximum number of spin 
flips. The summation of the corresponding ladder dia- 
grams can be performed by solving the integral equation 
for the c-f T-matrix (Fig. 0a)), 

T /<?%<?> ( iiJ n-,iu' n ,iCl) = -V 2 G b (iuj n + iJ n - itt)8 sc ,>5 s , a 
+ V 2 T^G b (iuj n + ioj[[ - in) x 

G /CT (i<)G c 4-< + ffl)T^5 v ,(<,^;,in). (4) 

Inserting NCA Green's functions for the intermediate 
state propagators of Eq. (||), we find numerically at low 
temperatures a pole of T^ c ^ in the singlet channel as a 
function of the center-of-mass (COM) frequency in the 
Kondo regime (rid > 0.7) pl| , p7t . This signals the ten- 
dency to form a singlet bound state at f2 = tt c f ~ — Tk- 
In the empty orbital regime (nd —> 0) the behavior of 
the system is governed by charge fluctuations. The dom- 
inating contributions in this low density region may be 
expected to result from conduction electron boson scat- 
tering. The corresponding scattering amplitude is 
obtained from Eq. (Q) by interchanging pseudofermions 
and antibosons, again leading to a pole, at £l cb < 0. In 
the mixed valence regime (rid — 0.5), the poles in both 
T(^) and T( cfc ) are of equal importance. 

In order to guarantee gauge invariance, self- 
consistency has to be imposed. The self-energies £/ , £& 
calculated from T^ c ^ and T' cfc ) then follow from a gen- 
erating functional <& [|l7| and are depicted in Fig. |l]b),c). 
They are given as nonlinear and nonlocal (in time) func- 
tionals of the Green's functions. The Green's functions 
in turn are expressed in terms of the self energies, closing 
the set of self-consistent equations (conserving T-matrix 
approximation, CTMA). Note that the contribution to $ 
containing one boson rung corresponds to NCA. The dia- 
gram with two rungs is excluded since it is not a skeleton. 
The sum of the $ diagrams with up to four rungs consti- 
tutes a large N expansion correct up to 0(lJN 2 ) and is 
identical to the diagram class used in Ref. 122] . We em- 
phasize that the CTMA, i.e. the self consistent summation 
of the infinite series of all diagrams shown in Fig. [I] is jus- 
tified on physical as well as formal grounds: At any loop 
order of $ it includes (1) the maximum number of spin 
flip as well as charge fluctuation processes; (2) all leading 
and sub-leading IR singular contributions, because all 
terms not included cancel pairwise in the IR regime [p3| . 
The threshold property of the auxiliary spectral func- 
tions implies that the exact T-matrices T^ c ^ and 
have no spectral weight at negative COM frequencies O, 
in contrast to the poles appearing in the "perturbative" 
evaluation, i.e. inserting NCA propagators as discussed 
after Eq. (||). Consequently, these poles are shifted to 
SI = by self-consistency, where they merge with the 
continuous spectral weight present for Q > 0, thus renor- 
malizing the threshold exponents of the auxiliary 
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FIG. 2. Pseudofermion and slave boson spectral functions 
Af and Ab in the Kondo regime (N — 2; Ed. = —0.05, 
r = 0.01 in units of the half-bandwidth), for a) the sin- 
gle-channel (M = 1) and b) the multi-channel (M = 2) case. 
In a) the symbols represent the results of NRG for the same 
parameter set, T = 0. The slopes of the dashed lines indi- 
cate the exact threshold exponents as given by Eqs. (2), (3) 
for M = 1 and by CFT for M = 2. Deviations from the power 
laws at low frequencies u) shift towards u) = as T — > 0, i.e. are 
finite T effects. Insets show magnified power law regions. 



spectral functions, as seen below. This is an expression 
of the fact that the Kondo singlet is not a two-particle 
bound state but rather a collective many-particle state. 

After analytical continuation to the real frequency axis 
we have solved the CTMA numerically by iteration. In 
the Kondo regime (n d > 0.7) of the N = 2, M = 1 
model, we have obtained reliable results down to tem- 
peratures of the order of 10~ 2 Tk (Note that Tk — » 
in the Kondo limit). In the mixed valence and empty 
impurity regimes, significantly lower temperatures may 
be reached, compared to the low temperature scale of 
the model. As is shown in Fig. |]a), the spectral func- 
tions obtained are in good agreement with the results of 
NRG (zero temperature results), given the uncertainties 
in the NRG at higher frequencies. Typical behavior in the 
Kondo regime (Fig. ||a)) is recovered: a broadened peak 
in Ab at lu ~ \Ed\, representing the hybridizing d-level 
and a structure in Af at uj ~ Tk ■ Both functions display 
power law behavior at frequencies below Tk, which at 
finite T is cut off at the scale u> ~ T. The exponents 
extracted from the frequency range T < uj < Tk of our 
finite T results compare well with the exact result also 
shown (see insets of Fig. ||) . A similar analysis has been 
performed for a number of parameter sets spanning the 



3 









" a b NCA 






- « f NCA c 













0.0 0.2 0.4 0.6 0.8 1.0 

rid 

FIG. 3. CTMA results (symbols with error bars) for the 
threshold exponents a/ and at of Af and Ab, N = 2, M = 1. 
Solid lines: exact (Eqs. (2), (3)), dashed lines: NCA results. 



consistency. There cannot be a renormalization of the 
exponents in any hnite order self-consistent summation. 
A numerical evaluation of the CTMA yields good agree- 
ment with the known exact values in the single-channel 
case, indicating that CTMA recovers the FL behavior. 
By contrast, in the multi-channel case the singular con- 
tributions are ineffective in the limit — > 1 of the two- 
channel Anderson model, i.e. the non-FL state persists 
and the exponents known from CFT are recovered. 
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part by DFG. Part of the work by J.K. was performed at 
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complete range of d-level occupation numbers rid- The 
extracted power law exponents are shown in Fig. to- 
gether with error bars estimated from the finite frequency 
ranges over which the fit was made. The comparatively 
large error bars in the mixed valence regime arise because 
here spin flip and charge fluctuation processes are of equal 
importance, inhibiting the convergence of the numerical 
procedure. In this light, the agreement with the exact 
results is very good. This is evidence that the CTMA re- 
covers the signature of FL behavior present in the exact 
auxiliary particle dynamics of the single-channel Ander- 
son model. As a consequence of the conserving scheme, 
the FL behavior should be reproduced in physical prop- 
erties as well, when derived from the same generating 
functional These evaluations are in progress. 

In the multi-channel case with M > 1, N = M, 
it follows from the symmetry of the model under the 
transformation f a — > b\, b m — > Ed — ► —Ed that 
ctf,b(nd) = «b,/(l — rid)- For nd = 1, a/,6 are known 
from CFT (see above). It follows that NCA yields the 
exact exponents both in the Kondo and in the empty 
impurity limits of the multichannel model. At present, 
it is not known whether this is the case for arbitrary nd- 
Using NCA Green's functions as discussed above, we find 
again a pole in the spin singlet channel of the T-matrix 
of pseudofermions (/ s , s =f, J.) and conduction electrons 
of flavor to (c CTm ) for nd < 1. However the weight of 
the pole vanishes for nd — ► 1, and a numerical solution 
of the self-consistent CTMA equations for M = N = 2, 
n d = 0.877 indeed yields exponents af ~ 0.44, a b ~ 0.49, 
very close to the exact value of 1/2 for rid = 1 (Fig. ||b))- 
This is consistent with the fact that there should not be 
a bound state contribution in the overscreened case. 

In this Letter we have considered the dynamics of the 
auxiliary particles for the (N,M) generalized Anderson 
impurity model with particular emphasis on the realiza- 
tion of FL behavior. We have shown that the occurrence 
of FL behavior can be deduced from the IR threshold 
exponents of the auxiliary particle spectral functions. 
A conserving self-consistent approximation incorporat- 
ing an infinite number of coherent spin flip and charge 
transfer processes (CTMA) leads to singular contribu- 
tions which renormalize the threshold exponents by self- 
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